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A one dimensional, parity-time ("PT)-symmetric magnetic metamaterial comprising split-ring 
resonators having both gain and loss is investigated. In the linear regime, the transition from the 
exact to the broken "PT-phase is determined through the calculation of the eigenfrequency spectrum 
for two different configurations; the one with equidistant split-rings and the other with the split- 
rings forming a binary pattern (VT dimer chain). The latter system features a two-band, gapped 
spectrum with its shape determined by the gain/loss coefficient as well as the inter-element coupling. 
In the presense of nonlinearity, the VT dimer chain with balanced gain and loss supports nonlinear 
localized modes in the form of novel discrete breathers below the lower branch of the linear spectrum. 
These breathers, that can be excited from a weak applied magnetic field by frequency chirping, can 
be subsequently driven solely by the gain for very long times. The effect of a small imbalance 
between gain and loss is also considered. Fundamendal gain-driven breathers occupy both sites 
of a dimer, while their energy is almost equally partitioned between the two split-rings, the one 
with gain and the other with loss. We also introduce a model equation for the investigation of 
classical VT symmetry in zero dimensions, realized by a simple harmonic oscillator with matched 
time-dependent gain and loss that exhibits a transition from oscillatory to diverging motion. This 
behavior is similar to a transition from the exact to the broken VT phase in higher-dimensional 
VT— symmetric systems. A stability condition relating the parameters of the problem is obtained in 
the case of piecewise constant gain/loss function that allows for the construction of a phase diagram 
with alternating stable and unstable regions. 
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PACS numbers: 



I. INTRODUCTION 



The investigation of artificial materials whose proper- 
ties can be tailored has recently attracted a lot of at- 
tention. Considerable research effort has been invested 
in the developement of artificial structures that exhibit 
properties not found in nature. Two recent and well 
known paradigms are the metamaterials, that provide 
full access to all four quadrants of the real permittiv- 
ity - permeability plano^, and the parity - time (VT) 
symmetric systems, whose properties rely on a delicate 
balance between gain and loss. The latter belong to a 
class of 'synthetic' materials that do not obey separately 
the parity (V) and time (T) symmetries but instead 
they do exhibit a combined VT symmetry. The ideas 
and notions of VT— symmetric systems have their roots 
in quantum mechanics where VT— symmetric Hamilto- 
nians have been studied for many yearsi^. The notion 
of VT symmetry has been recently extended to dynam- 
ical lattices, particularly in optics, where photonic lat- 
tices combining gain and loss elements offer new possi- 
bilities for shaping optical beams and pulses. Soon af- 
ter the developement of the theory of VT— symmetric 
optical lattices^^, the VT— symmetry breaking was ex- 
perimentally observe d 12 ' 29 ' 36 . Naturally, such consider- 
ations have been also extended to nonlinear Iattices2i24 
and oligomers 2 ^, and VT— related phenomena like unidi- 



rectional optical transport^!, unidirectional invisibility^ 2 -, 
and Talbot effects 2 ^ were theoretically demonstrated. 
Moreover, it has been shown that optical solitonsi 2 -^^, 
nonlinear modes^ 9 -, and breather a 4 ' 20 may also be sup- 
ported by VT— symmetric systems. Moreover, the ap- 
plication of these ideas in electronic circuits^, not only 
provides a platform for testing VT— related ideas within 
the framework of easily accessible experimental configu- 
rations, but also provides a direct link to metamaterials 
whose elements can be modeled with equivalent electrical 
circuits. 

Conventional metamaterials comprising resonant 
metallic elements operate close to their resonance 
frequency where unfortunately the losses are untolerably 
high and hamper any possibility for their use in device 
applications. The pathways to overcome losses are either 
to replace the metallic parts with superconducting ones 
(superconducting metamaterials)^, or to construct active 
metamaterials by incorporating active constituents 
that provide gain through external energy sourses. 
The latter has been recently recognized as a very 
promising technique for compensating losse o 5 ' 6 ' 33 . A 
particular electronic component that may provide both 
gain and nonlinearity in a metamaterial is the tunnel 
(Esaki) diode which features a current-voltage charac- 
teristic with a negative resistance part£. Left-handed 
transmission lines with successful implementation of 



Esaki diodes have been recently realized^, although 
other electronic components may be employed as 
well for loss compensation^. Thus, the fabrication 
of VT— symmetric metamaterials with balanced gain 
and loss is feasible with the present technology in the 
microwaves, combining highly conducting split-ring 
resonators (SRRs) and negative resistance devices in 
a way similar to that in electronic circuits^. In this 
prospect, the SRR equivalent circuit parameters and the 
bias of the negative resistance device sould be properly 
adjusted to provide gain and equal amount of loss, as 
well as real eigenfrequencies in a finite frequency range 
of the gain/loss parameter. 

In the following we present a one-dimensional, discrete, 
equivalent circuit model for an array of SRRs with alter- 
natingly gain and loss in the two different configurations 
(Section II). In Section III we present linear eigenfre- 
quency spectra for systems with small number of SRRs 
and we obtain the relation that provides the eigefrequen- 
cies for large systems. It is shown that VT— symmetric 
metamaterials undergo spontaneous symmetry breaking 
from the exact VT phase (real eigenfrequencies) to the 
broken VT phase (at least a pair of complex eigenfre- 
quencies), with variation of the gain/loss coefficient. In 
Section IV, where nonlinearity becomes important, the 
generation of long-lived nonlinear excitations in the form 
of discrete breathers (DBs)^i is demonstrated numeri- 
cally. These novel gain-driven DBs result by a purely 
dynamical proccess, through the matching of the input 
power through the gain mechanism and internal loss. In 
Section V we introduce a model VT— symmetric system 
in zero dimensions, realized by a harmonic oscillator with 
balanced time-periodic gain and loss, that exhibits ex- 
traordinary properties and multiple critical (phase tran- 
sition) points. Section VI contains the conclusions. 
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FIG. 1: Electrical equivalent circuit for a split-ring resonator 
loaded with a tunnel (Esaki) diode. 



II. EQUIVALENT CIRCUIT MODELLING AND 
DYNAMIC EQUATIONS 

Consider a metallic split-ring resonator (SRR) that can 
be regarded as an RLC electrical circuit featuring Ohmic 
resistance R, inductance L, and capacitance C. A tun- 
nel (Esaki) diode is connected in parallel with the ca- 
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FIG. 2: (Color online) Schematic of a VT— symmetric 
metadimer comprising two tunnel diode-loaded SRRs in an 
alternating magnetic field H(t). The SRRs are coupled mag- 
netically through their mutual inductance M. Different bias 
in the diodes may create a balanced gain/loss structure. 



pacitance C of the SRR (Fig. 1) forming thus a non- 
linear metamaterial element with gain. Esaki diodes ex- 
hibit a well defined negative resistance region in their 
current-voltage characteristics that has a characteristic 
'N' shape. A bias voltage applied to the diode can move 
its operation point in the negative resistance region and 
then the SRR-diode system gains energy from the source. 
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FIG. 3: (Color online) Schematic of a one-dimensional 
VT— symmetric metamaterial. Upper panel: all the split- 
ring resonators are equidistant. Lower panel: the separation 
between split-ring resonators is modulated according to a bi- 
nary pattern {VT dimer chain). The applied field is such that 
its magnetic component is perpendicular to the plane of the 
split-rings. 

A metadimer comprising two SRRs loaded with tun- 
nel diodes in an external alternating magnetic field is 
shown in Fig. 2. The equivalent circuit parameters 
R, C, and L of the SRRs and the bias in the diodes 
have been adjusted so that: (i) the two elements have 
the same eigenfrequencies; (ii) one of the SRRs has gain 
while the other has equal amount of loss. Then, the pair 
of SRRs is a VT— symmetric metadimer that can used 
for the construction of a one-dimensional VT— symmetric 
metamaterial, which moreover are nonlinear due to the 
tunnel diodes. The alternating magnetic field induces 
an electromotive force (emf) in each SRR due to Fara- 
day's law which in turn produce currents that couple the 



SRRs magnetically through their mutual inductance M 
(Fig. 2). The coupling strength between SRRs is rather 
weak due to the nature of their interaction (magnetoin- 
ductive), and has been calculated accurately by several 
author a 28 ' 35 . The SRRs may also be coupled electrically, 
through the electric dipoles that develop in their slits. 
Thus, in the general case one has to consider both mag- 
netic and electric coupling between SRRs. However, for 
particular relative orientations of the SRR slits the mag- 
netic interaction is dominant, while the electric interac- 
tion can be neglected in a first approximatio n 10 ' 13 ' 32 . As 
can be seen in Fig. 3, the VT— symmetric metadimers 
can be arranged in a one-dimensional lattice in two dis- 
tinct configurations; one with all the SRRs equidistant 
and the other with the SRRs forming a VT dimer chain. 
Within the framework of the equivalent circuit 
model, a set of discrete differential equations has been 
used to describe the dynamics in nonlinear magnetic 
metamaterial a 16 ' 19 ' 25 ' 28 . Takining into account the bi- 
nary structure of the VT dimer chain, the dynamics 
of the VT— symmetric metamaterial with balanced gain 
and loss is governed by the following equations which are 
presented in normalized form^S 



^M?2n + <?2n+l + Am<727i+2 
+ X E q2n + <?2n+l + ^E<iln+1 

- a( lln+i + P<lln+i + l<J2n+i = £o srn(fir) 

Am<Z2ti-1 + Q2n + A M <72n+l 
FAg^n-l + l2n + A s <72n+1 

792n = £o sin(fir) 
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where Xm , X' M and Xe , X' E are the magnetic and elec- 
tric coupling coefficients, respectively, with Xe.m > X' E M 
and A^mA^ m > 0, a and /3 are dimensionless nonlin- 
ear coefficients, 7 is the gain/loss coefficient (7 > 0), £0 
is the amplitude of the external driving voltage, while 
fl and r are the driving frequency and temporal vari- 
able, respectively, normalized to the inductive-capacitive 
(LC) resonance frequency ujq and inverse LC resonance 
frequency 0Jq , respectively, uio = 1/y/LCa with Co being 
the linear capacitance. The values selected for the non- 
linear coefficients a = —0.4, f3 = 0.08 are typical for a 
diode and they provide a soft on-site nonlinear potential 
for each metamaterial element. They can be obtained 
from a Taylor expansion of the capacitance-to -oltage re- 
lation of an equivalent circuit diode model, that gives a 
very good approximation for weakly driven system o 19 ' 37 . 



III. LINEAR EIGENFREQUENCY SPECTRA 
AND CRITICAL POINT 

In order to obtain the critical value of 7 = 7 C that sep- 
arates the exact VT— phase, where all the eigenvalues are 
real, from the broken VT— phase, where at least a pair of 
eigenvalues is complex, we calculate the frequency spec- 
trum. This is a straightforward procedure for systems 




FIG. 4: (Color online) Frequency eigenvalues of the free 
VT— symmetric SRR array as a function of the gain/loss pa- 
rameter 7 for X E = 0, X M = -0.1, and (a) N = 2; (b) JV = 4; 
(c) N = 8; (d) N = 16. The arrows indicate the critical point 
7c for the exact-to-broken VT phase transition. Only the real 
eigenfrequencies is shown for clarity. 



with relatively small number of SRRs; the roots of the 
determinant of the linearized Eqs. (JlJ and ((2|) for £0 = 
are obtained with a root-finding algorithm and then plot- 
ted against the gain/loss parameter 7. In Figs. 4 and 5, 
the real eigenfrequencies of VT— symmetric metamate- 
rials in both configurations are shown as a function of 
7, while the arrows indicate the critical point j c in each 
case. Thus, for 7 < -f c all eigenvalues are real, while in 
the opposite case at least a pair of eigenvalues has become 
complex. As we can see from the figures, for 7 > 7 C more 
and more eigenfrequency pairs become complex with in- 
creasing 7, until they all become complex for a particular 
value of 7. Moreover, as we can see from an inspection of 
Figs. 4 and 5, obtained for the equidistant SRR configu- 
ration and the VT dimer chain, respectively, the value of 
7 C decreases rapidly with increasing number of SRRs ./V 
for the former configuration, while it tends to a constant 
value for the latter configuration. This can be seen more 
clearly in Fig. 6, where the critical point 7 C is plotted as 
a function of N for both configurations. For the curves 
corresponding to equidistant SRRs (corresponding to two 
different values of the magnetic coupling coefficient Am) 
we see that 7 C is smaller for lower magnetic coupling co- 




FIG. 5: (Color online) Frequency eigenvalues of the free 
VT— symmetric dimer chain as a function of the gain/loss 
parameter 7 for \e — A' E = 0, Am — —0.1, Am = —0.05, 
and (a) N = 2; (b) N = 4; (c) N = 8; (d) N = 16. The ar- 
rows indicate the critical point y c for the exact-to-broken VT 
phase transition. Only the real eigenfrequencies are shown for 
clarity. 



efficient Am- However, in both curves corresponding to 
equidistant SRRs the value of the critical point 7 C tends 
to zero with increasing N. In contrast, for the VT dimer 
chain configuration, the value of 7 C tends to a constant 
finite value which approximatelly equals the absolute dif- 
ference of the magnetic coupling coefficients Am and A' M 
(see below). 

For large systems, we can obtain a condition that de- 
termines the critical point "f c as a function of the mag- 
netic coupling constant(s). In the standard way, we sub- 
stitute into the linearized Eqs. (JTJ and (J2J for Eq = the 
trial solutions 

q 2n = Aexp[i(2nK — Or)], (3) 

q 2n +i = Bcxp{i[(2n + l) K ~nT}}, (4) 

where k is the normalized wavevector. Then, by request- 
ing nontrivial solutions for the resulting stationary prob- 
lem, we obtain 



ft 2 = (-b±y/b 2 -4ac\ /(2a), 



(5) 




FIG. 6: (Color online) Dependence of the critical gain/loss 
parameter value y c on the number of SRRs, N, for mag- 
netically coupled SRRs in both the equidistant and dimer 
chain configuration. The black squares and the red dia- 
monds have been calculated for the former configuration with 
Am = —0.1 and Am — —0.05, respectively. The green cir- 
cles have been calculated for the latter configuration with 
Am = —0.1, A'm = —0.05. The lines serve as a guide to the 
eye. 



where 



0= 1 
-2[1- 



- (A M 
(A b - 



V, 



Mi 



Ae)(Am - A' M )] 
c = I - (A E - A' E ) 2 -s K e' K , 
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2Aecos(k), 



2A^cos(k), /j, k = 2Amcos(k), 



and t 

fj,' K = 2A^ 7 cos(k). In the following, we neglect the electric 



coupling between SRRs, i.e., Ae 
Then, Eq. ([5]) reduces to 



A^ = 0, for simplicity. 



nl = 



2 - 7 2 ± ^ 7 4 ~ 2 7 2 + (Am - A^) 2 + ^g 
2(1 - (Am - X' M ) 2 - MkM'k) 



(9) 



The condition for having real fl K for any k then reads 



cos 2 (k) > 



7 2 (2- 7 2 )-(A 
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It is easy to see that for Am 
to the equidistant SRR configuration, the earlier con- 
dition cannot be satisfied for all k's for any positive 
value of the gain/loss coefficient 7, implying that a large 
VT— symmetric SRR array (Fig. 3, upper panel) will be 
in the broken phase. To the contrary, for Am 7^ A' M , 
i.e., for a VT dimer chain (Fig. 3, lower panel), the 
above condition is satisfied for all re's for 7 < 7 C — 
I Am — A M |, (7 4 ~ 0). In the exact phase (7 < j c ), the 
VT— symmetric dimer array has a gapped spectrum with 
two frequency bands, as shown in Fig. 7. The width of 
the gap separating the bands decreases with decreasing 
I Am — A M | for fixed 7. For 7 ~ j c the gap closes, some 
frequencies in the spectrum acquire an imaginary part 
and the VT metamaterial enters into the broken phase. 
Note that the gain/loss coefficient 7 has little effect on 
the dispersion curves of the VT dimer chain (compare 




FIG. 7: (Color online) Frequency bands for a VT— symmetric 
dimer chain with balanced gain and loss for Am = —0.17, 
\' M = —0.10, and 7 = 0.05 (black solid curves); 7 = (red 
dotted curves). The gaps are indicated in green (dark) color. 



with the dotted curves where 7 is set to zero), as long as 
the sign in front of 7 alternates from one SRR to another. 



IV. GAIN-DRIVEN BREATHER EXCITATIONS 

For a gapped linear spectrum, large amplitude linear 
modes become unstable in the presence of driving and 
nonlincarity. If the curvature of the dispersion curve in 
the region of such a mode is positive and the lattice po- 
tential is soft, large amplitude modes become unstable 
with respect to formation of DBs in the gap below the 
linear spectrum^. For the parameters used in Fig. 7, 
the bottom of the lower band is located at fio — 0.887, 
where the curvature is positive. The corresponding pe- 
riod at the bottom of the lower band is To = 27r/Oo- 
Moreover, the SRRs are subjected to soft on-site poten- 
tials for the selected values of the nonlinear coefficients 
a and f3. Then, DBs can be generated spontaneously 
by a frequency chirped alternating driver; after turning 
off the driver, the breathers are driven solely by gain. 
A similar procedure has been applied succesfuly to lossy 
nonlinear metamaterials with a binary structur e) 17 ' 18 ' 25 
The results are illustrated in Figs. 8 and 9, where the 
case of a slight imbalance between gain and loss and its 
effect on breather generation has been also considered. 
In these figures, a density plot of the local energy E n of 
a VT— symmetric metamaterial is shown on the n — r 
plane for two different values of the driving amplitude 
e . 

In the following, the integration of Eqs. ([T]) and ([2]) 
implemented with the boundary condition 



Qo(r) = qN+i(r) = 0, 



(11) 



that accounts for the termination of the structure in fi- 
nite systems, is performed with a 4th order Runge-Kutta 
algorithm with fixed time-step. In order to prevent in- 
stabilities that will result in divergence of the energy at 



particular sites in finite time, we consider a longer dimcr 
chain with total number of SRRs N + 2Nf, then we re- 
place the gain with equivalent amount of loss at exactly 
Ng SRRs at each end of the extended chain. In other 
words, we embbed the VT— symmetric dimer chain into 
a geometrically identical lossy chain, in order to help the 
excess energy to go smothly away during evolution liv- 
ing behind stable (or at least very long-lived) breather 
structures. 

We use the following procedure described in detail in 
Ref.^: 

• At time r = 0, we start integrating Eqs. ([1]) and 
@ from zero initial state without external driving for 
500 Tq ~ 3500 time units (t.u.) to allow for significant 
developement of large amplitude modes. 

• At time t ~ 3500 t.u. the driver is switched-on 
with low-amplitude Eq and frequency slightly above Qo 
(1.01 Slo — 0.894). The frequency is then chirped down- 
wards with time to induce instability for the next 10600 
t.u. (~ 1500 T ), until it is well below fi (0.997 £7 ~ 
0.882). During that phase, a large number of excitations 
are generated that move and strongly interact to each 
other, eventually merging into a small number of high 
amplitude breathers and multi-breathers. 

• At time r ~ 14100 t.u. (point A on Figs. 8 and 9), the 
driver is switched off and the DBs that have formed are 
solely driven by the gain. They continue to interact for 
some time until they reach an apparently stationary state 
and get trapped at particular sites. The high density 
segments between the points A and C in Figs. 8 and 9 
precisely depict those gain-driven DB structures. 

• At time r ~ 440000 t.u. (point C on Figs. 8 and 9), 
the gain is replaced everywhere by equal amount of loss, 
and the breathers die out rapidly. Note that the above 
procedure of breather generation is very sensitive to pa- 
rameter variations of the external fields. Even though 
the values of the driving amplitudes in Figs. 8 and 9 are 
rather close (i.e., £0 = 0.085 and 0.095, respectively), the 
breather structures as well as their numbers are differ- 
ent. In Fig. 8(b), in the balanced gain/loss case, we ob- 
serve two distinct structures that have been formed that 
correspond to a relatively high amplitude multi-breather 
and a low amplitude breather. These structures remain 
stationary during the long time interval they have been 
followed (> 56000 T ). In Figs. 8(a) and 8(c), gain and 
loss are not perfectly matched; in Fig. 8(a) loss exceeds 
gain by a small amount while in in Fig. 8(c) gain exceeds 
loss by the same small amount. Notably, breather excita- 
tions may still be formed through the frequency chirping 
procedure in the presence of a small amount of cither net 
gain or net loss. Indeed, as we may observe comparing 
Figs. 8(a) and 8(c) with Fig. 8(b), the same structures 
are formed [except the low amplitude breather that is not 
visible in Fig. 8(c)]. However, both in Figs. 8(a) and 8(c) 
we can see a slow gradual widening of the high amplitude 
multibrcathcr: when loss exceeds gain the multibrcathcr 
losses its energy at a low rate, with its excited sites that 
are closer to its end-points gradually falling down to a 





FIG. 8: (Color online) Spatiotemporal evolution of the energy 
density E n for a VT— symmetric dimer chain with N = 70, 
N e = 10, Q = 0.887, 7 = 0.002, Am = -0.17, X' M = -0.10 
(Xe — X' E — 0), so = 0.085, 7 = 0.002 and (a) excess loss 
0.2%; (b) balanced case; (c) excess gain 0.2%. 



FIG. 9: (Color online) Spatiotemporal evolution of the energy 
density E n for a VT— symmetric dimer chain with N = 70, 
N t = 10, Qo = 0.887, 7 = 0.002, Am = -0.17, X' M = -0.10 
(Xe — X' E — 0), eo = 0.095, 7 = 0.002 and (a) excess loss 
0.2%; (b) balanced case; (c) excess gain 0.2%. 



low amplitude state. Similarly, when gain exceeds loss, 
the high amplitude multibreather slowly gains energy and 
becomes wider. In both cases, breather destruction will 
take place in a time-scale that depends exponentially on 
the gain/loss imbalance. Thus, in an experimental situa- 
tion, where gain/loss balance is only approximate, it will 
be still possible to observe breathers at relatively short 
time-scales. Similar observations hold for Fig. 9 as well. 
In this figure, we observe three relatively high amplitude 
multibreather structures that are formed both in the bal- 
anced and the imbalanced case. Here we also observe 
that instabilities may appear after long time intervals of 
apparently stationary breather evolution as well. When- 
ever this happens, they start moving through the lat- 



tice until they get once more trapped at different lattice 
sites. In Fig. 9(b) such an instability appears between 
200000 — 250000 t.u. for the two narrower multibreather 
structures. The one of them gets trapped a few tenths 
lattice sites away from its previous position, while the 
other (the narrowest) collides and it is absorbed by the 
wide multibreather located at n 50. 



V. CLASSICAL VT -SYMMETRY IN ZERO 
DIMENSIONS 

In the classical domain, in all cases of VT— symmetric 
systems investigated so far, the combination of time- 



reversal and parity symmetries is utilized: In a given 
time-evolution if we reverse time while "reflecting" posi- 
tion and momentum, we must retrace the original path. 
The parity operation requires that the system is extended 
either in continuous or discrete space in order to be able 
to perform this operation. We show here that this is not 
necessary and that the features of the VT -symmetric 
systems can be preserved in "zero" dimensions where the 
parity symmetry is trivial. 

We consider the following simple harmonic oscillator: 



29{t)i 



uj x = 0, 



(12) 



where x = x(t) is the equilibrium dislacement of a mass 
(the charge in an RLC cirquit), ujq the resonant oscilla- 
tion frequency while 9(t) is a time-dependent "damping" 
term; both ljo and 9(t) are scaled to the mass (impedance 
of the circuit). We take 9(t) to be periodic with period T, 
viz. 9(t + T) = 9 {t) while its values may be both positive 
and negative, i.e., for some part of the period the oscil- 
lator experiences damping while in the rest of the period 
time amplification, or anti-damping. We investigate the 
osillator evolution after long time and the stability of the 
motion. Instead of addressing a general periodic function 
9(t) we focus on a simple form that makes the problem 
readily solvable, viz. we take 



9{t) 



+7 if < t < n; 

— 7 if Ti < t < T%. 



(13) 



where T = t\ + t 2 and taking the plus (minus) sign in 
front of the coefficient 7 (7 is defined as a positive con- 
stant, 7 > 0, that may assume any value between zero 
and unity) we have in the first (second) part of the cycle 
loss (gain). With this form of piecewise constant function 
9(t) we can easily solve Eq. (fT2f for the loss (L) segment 
of time duration T\ and gain (G) segment of duration ti- 
The form of Eq. (|13p permits to view the problem as 
mapping of the position- velocity vector at a given time 
to the position-velocity vector at a later time; if in the 
begining of the gain (loss) period (assuming t = 0) we 
have position and velocity equal to (.To, Xq). Then after 
the evolution during time r (n or t 2 , respectively) we 
obtain: 



Mc/dr) ( 2 



(14) 



where for gain we have 



olT! 



M G (n) = — 



8 cos 8t\ — 7 sin 8t\ sin 8t\ 

— luq sin 5t\ 5 cos 5t\ + 7 sin 8t\ 

and for the loss respectively 



(15) 



_ 7T2 



M L (n) : 

6 cos 8t 2 + 7 sin 8t 2 sin 8r 2 

—cuq sin 8r 2 8 cos 8t\ — 7 sin 8t 2 



(16) 




FIG. 10: (Color online) Phase diagram on the 7 — r plane for 
the harmonic oscillator with time-dependent " damping" term, 
obtained from Eq. (|27|l . The stability region is indicated in 
blue (dark) color. 



where 8 = \fu 



7 2 . Using the mapping, we may ob- 



tain long time evolution after N periods T as a repeti- 
tive application of the matrices Mq(t{) and Ml(t 2 ) to 
an arbitraty initial state (xo, xq). Since the matrices 
A/L/G( T ) e±7T are unimodular, the long time evolution 
with be dominated by the exponential term exp[7V7(Ti — 
t 2 )]\ this leads trivially to exponential growth (n > r 2 ) 
or exponential decay (t\ < t 2 ). As a result we consider 
the more interesting case with t\ = t 2 = r = T/2; in 
this case the gain and loss power is perfectly matched 
during the period T. The combined propagation ma- 
trix after one period (assuming first gain) is simply 
M{T) = M l (t)M g (t), i.e.,' 



M(T) 



J_ / Mn M12 
£2 I M 2 i M 22 



where Mj ? - are given by 



M 
M 21 



Mn = -7 2 +WqCos(20), 

= +2 sin 4>(8 cos + 7 sin </>) 

-2w sin <f>{8 cos <fi — 7 sin cj)) 

M 22 = -7 2 + ul cos(20) 



(17) 



(18) 
(19) 
(20) 
(21) 



and 4> = St = ST/2. The matrix M is clearly also uni- 
modular with eigenvalues e v and e~ % ^\ since the trace of 
a matrix is invariant we find 



cos/i = 
Stability Equation.- Eq. 



l-B 2 cos2< ? (> 
cos ^ = — 7^2 = 1 



-7 2 + loq cos 2(j> 

can be re written as 
2B 2 sin 2 6 



l-B 2 



(22) 



(23) 



with B = — . In order to have stable motion it is neces- 
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sary that | cos/i| < 1, or — 1 < 1 -\ 1 _ S g2 > 1, leading 



to 



„ 2B 2 sin 2 (ft 

2 < < 1 

" 1 - B 2 ~ 



(24) 



Eq. ([24]) has solutions only for \B\ > 1 or |— ^ | < ojq, for 
loq > 0; in the latter case we find 



> 



UlQ 



(25) 



The range thus of allowed values for the angle (ft (7 > 0) 

is 



-1 < COS 



< COS (ft < 1 



< 



_1_ 

Wo 



(26) 



We note that there are multiple allowed solutions marked 
by the the lines cos0 = ±7/^0, i.e., for the three param- 
eters of the problem 7, T = 2r and ujq the transcedcntal 
equation 



cos((5r) 



+3- for 2mr - tt/2 < St < 2mr + tt/2: 
—3- for 2nn + tt/2 < St < 2tt?i + 3tt/2, 



where n — 0,±1,±2, .., marks the onset of the transition 
from stable to unstable evolution. This is the equivalent 
to the VT transition from the exact (stable) to the bro- 
ken (unstable) phase in this zero-dimensional problem. 
Using Eq. (f25j) we construct a VT phase diagram on the 
r — 7 plane (Fig. 10), where the blue (dark) color indi- 
cates regions of stability. If we fix one of the parameters, 
variation of the other drives the oscillator through alter- 
natingly stable and unstable regions, as can be readily 
observed from Fig. 10. 




FIG. 11: (Color online) Density plot of the logarithm of the 
frequency spectra of x(t), y = log 10 {PS[x(t)]}, as a function 
of 7. The discrete frequency components of the spectra for 
each value of 7 are indicated with blue (darker) color. 

Introducing the frequency u = 2tt/T = tt/t we define 
the reduced parameters tt = lu/ujq and 7 = j/ujq. Then, 



the equation that controls the stability regions becomes 

(27) 



cos 



TiV 1 ~7 2 



O 



1, 



where the tilde has been dropped. Eq. (|2"T|) has solutions 
for I7I < 1; for different reduced frequencies £1 we obtain 
different number of solutions of the earlier equation, the 
number of the latter icreases with decreasing frequency 
co. Once we have the solutions of Eq. (|2"7| we can find the 
regimes of stability and instability. Consider the resonant 
case f2 = 1 where the external frequency of gain/loss 
alternation matches the self-frequency of the oscillator. 
Besides the trivial solution at 7 = 0, numerical solution 
of Eq. (|2"T|) gives 71 w 0.676 so that the stable region is 
in the range < 7 < 7i- For fi = 0.5 the corresponding 
solutions are 71 ps 0.54, 72 w 0.80, and 73 w 0.90 with 
two stability regions, i.e., < 7 < 7i and 72 < 7 < 73. 

The stable solutions of Eq. (fT2")) as a function of time 
t are in general quasi-periodic oscillations whose spectral 
content varies both with r and 7. For an illustration, we 
consider a particular value of r, i.e., r = 2-k (fi — 0.5) for 
which the boundaries of the stability regions have been 
calculated. For the interval with relatively low values of 
7> < 7 < 7i — 0.54, that is more physically relevant, we 
present a density plot of the logarithm of the frequency 
spectra of x(t), y = log 10 {PS[x(t)]}, as a function of 7 
(Fig. 11). In this figure, the discrete frequency compo- 
nents of the spectra for each value of 7 are indicated in 
dark (blue) color. For 7 = and very close to zero, the 
only frequency appearing in the spectrum is the eigen- 
frcquency of the oscillator $7 ~ 1. With increasing 7, the 
frequency components at £1 ~ 1 ± 0.5 become more and 
more important. At about 7 = 0.15 these half-integer 
frequency components start splitting into pairs of fre- 
quencies that are symmetric around the half-integer val- 
ues. The separation between these pairs increases with 
further increasing 7, and the frequency components from 
neighboring pairs come closer and closer together until 
they eventually merge for 7 approaching its critical value 
0.54. 



VI. CONCLUSIONS 

We have investigated theoretically a VT— symmetric 
nonlinear metamaterial relying on gain and loss. Eigcn- 
frequency spectra for linearized systems of either small 
and large N and two different configurations were cal- 
culated and the critical points j c were determined. 
Large VT— symmetric metamaterials with the dimcr 
chain configuration exhibit phase transitions for the ex- 
act to the broken VT— symmetry phase, while large 
VT— symmetric metamaterials with the equidistant SRR 
configuration are always in the broken VT phase. 

In the presence of nonlinearity, we have demonstrated 
numerically the a VT— symmetric dimer chain supports 
localized excitations in the form of discrete breathers. 



Breathers are excited by a purely dynamical process, 
with a frequency chirped external magnetic field that in- 
duces instability in a zero initial state. Subsequently, 
the nonlinearity focuses energy around points that have 
acquired high amplitude leading to the formation of lo- 
calized structures. The external field is then switched off 
and those localized structures are then solely driven by 
the gain, while the excess energy lives the system through 
its lossy parts at the ends, leading eventually to breather 
generation. 

Remarkably, slight imbalance between gain and loss 
does not destroy the breathers instantly; they can still 
be generated through the frequency chirping procedure 
and they can be regarded as stationary for relatively 
short time intervals. In the long term, the imbalanced 
breathers either gain constantly energy and diverge or 
lose constantly energy and vanish. The time-scale for 
the latter events depends exponentially on the amount 
of imbalance. 

We have introduced a "zero-dimensional" VT system 
that can be realized by a harmonic oscillator with a 
"damping" term that provides balanced gain and loss, 



through alternation of the sign of the damping coeffi- 
cient. We consider a piecewise linear gain/loss function 
and obtain a stability condition, i.e., a relation between 
the parameters of the problem. We thus obtain a "phase 
diagram" with parameter regions where oscillatory (sta- 
ble) motion and diverging (unstable) motion occur. A 
crossing of the stability boundary marks the onset of a 
transition from stable to unstable evolution that is equiv- 
alent to the VT transition from the exact to the broken 
phase in this zero-dimensional problem. 
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